Deterministic creation of stationary entangled states by dissipation 
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We propose a practical physical system for creation of a stationary entanglement by dissipation without 
employing the environment engineering techniques. The system proposed is composed of two perfectly distin- 
guishable atoms, through their significantly different transition frequencies, with only one atom addressed by 
an external laser field. We show that the arrangement would easily be realized in practice by trapping the atoms 
at the distance equal to the quarter-wavelength of a standing-wave laser field and locating one of the atoms at 
a node and the other at the successive antinode of the wave. The undesirable dipole-dipole interaction between 
the atoms, that could be large at this small distance, is adjusted to zero by a specific initial preparation of the 
atoms or by a specific polarization of the atomic dipole moments. Following this arrangement, we show that 
the dissipative relaxation can create a stationary entanglement on demand by tuning the Rabi frequency of the 
laser field to the difference between the atomic transition frequencies. The laser field dresses the atom and we 
identify that the entangled state occurs when the frequency of one of the Rabi sidebands of the driven atom tunes 
to frequency of the undriven atom. It is also found that this system behaves as a cascade open system where the 
fluorescence from the dressed atom drives the other atom with no feedback. 

PACS numbers: 03.65.Yz, 03.67.Bg, 42.50.Dv 



Engineering a long-lived entanglement between atoms or 
ions located at small distances and interacting with a vacuum 
reservoir is a challenging problem in the quantum informa- 
tion science fl-lH]. The problem has already been the sub- 
ject of several papers but attention was concentrated mainly 
on the effect of the direct interatomic dipole-dipole coupling 
and on the role of dark states for that entanglement iq 420ll . 
Although closely spaced atoms can occupy an area of small 
sizes, which has the advantage of controlling and manipulat- 
ing the dissipative relaxation in the system, it is practically 
difficult to keep the closely spaced atoms at fixed positions. At 
small distances the dipole-dipole interaction is strong causing 
the atoms to move or rapidly oscillate around their equilib- 
rium positions |21fl . Therefore, many schemes have been pro- 
posed for entangling atoms or atomic ensembles separate d by 
a large distance, much larger than an optical wavelength 122- 
l28tl . In this case, the direct dipole-dipole interaction between 
the atoms can be ignored, but at the cost of entanglement sta- 
bility due to difficulties in controlling the dissipative relax- 
ation of atoms spread over large regions. The relaxation is a 
source of decoherence that leads to irreversible loss of infor- 
mation encoded in the internal states of the atoms and thus is 
regarded as the main obstacle in the generation of persistent 
entanglement. In contrast, an idea of "quantum-state engi- 
neering" by dissipation have been invented, where the dissi- 
pative relaxation can be used to drive atoms to a long-lived 
entangled state JMS. 

In the following we consider a simple physical system 
where a stationary entanglement can be created by a dissi- 
pative exchange of photons between two separate systems 
coupled to a decohering environment. The proposed system 
is composed of two non-identical atoms which decay in a 
vacuum reservoir with no coherent couplings and no hidden 



initially populated entangled states present. The atoms are 
placed at a small but finite distance r\i and are driven by a 
standing-wave laser field propagating in the direction paral- 
lel to the inter-atomic axis and are also coupled to the rest of 
the multi-mode radiation field which is in an ordinary vacuum 
state. Each atom is modeled as a two-level system (qubit) with 
the ground state and the excited state |e»), (i = 1, 2) sep- 
arated by a transition frequency u>i. We assume that in the 
absence of the laser field the atoms are perfectly distinguish- 
able that the frequency difference Ao = u>\ — L02 is quite large 
to prevent photons emitted by each atom to induce changes in 
the other atom. In addition, we assume that the atomic dipole 
moments pi are not necessary equal but are parallel to each 
other, and arrange the external driving of the system such that 
only one atom is coupled to the laser field. This coupling 
could be achieved in practice by placing one of the atoms at 
a node and the other at the successive antinode of the stand- 
ing wave. With this arrangement, we shall assume that only 
atom 2 is addressed by the laser field. 

The dynamics of the system is described by the density op- 
erator p, which in the interaction picture satisfies the master 
equation 



p = --[H + H d ,p]+L(p), 



(1) 



where 



Ho = h(A Q + A L )S z i + HA L SI + -hn (S+ + Si) , (2) 

is the Hamiltonian of the atoms and the interaction of the 
atom 2 with the laser field, 



H d = nn 12 (s? S2 + h.c. 



(3) 
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is the coherent dipole-dipole interaction between the atoms 
with f2i2 = Kelli2, and the Liouvillian C(p) describing the 
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dissipative relaxation of the system is of the form 1351 3 

2 

C{ P ) = -J2 ^ ( I s ? > s Jp1 + IP 15 ?' S J] ) W 
with 7x2 = 721 = — ImZ7i2, and 

U 12 = 



2 ^((l-cos 2 ,) — 



+ (1-3 cos 2 r?) 



(fcri 2 ) 2 (kr 12 ) 3 



e lkri2 . (5) 



Here, S^(S~ ) is the raising (lowering) operator of the ith 
atom and Sf describes its energy, Oo is the resonant Rabi fre- 
quency of the laser field at the position of the driven atom, 
Al = uj2 — wl is the detuning of the transition frequency 
of the driven atom from the laser field frequency ujl, and 
2~fu = 2~fi is the relaxation rate of the ith atom. 

The parameters SI12 and 712, symmetrical with respect to 
the exchange of the atoms (1 H 2) even if the atoms are not 
identical, describe mutual interaction of the atoms through 
the coupling to the common reservoir. They depend on the 
distance ri 2 between the atoms and the orientation 77 of the 
atomic dipole moments in respect to the interatomic axis. The 
parameters are of importance only in systems of atoms sepa- 
rated by small or intermediate distances relative to the radia- 
tion wavelength, fcrx2 < 1, and vanish at fcri 2 ^> 1. However, 
they can also vanish, but not simultaneously at some finite dis- 
tances. We shall use this fact to eliminate the dipole-dipole 
interaction 0x2 leaving only the dissipative interaction 712. It 
is easily verified that 0x2 vanishes at the distance ri2 = A/4 
and for a fixed r\ = 0.3047T. In practical terms, this could 
correspond to atoms trapped at fixed positions, the Raman- 
Nath approximation, and the atomic dipole moments polar- 
ized at a fixed direction. This is compatible with many exper- 
iments on cooling of trapped atoms, where the storage time 
of the trapped atoms is long, so that they are essentially mo- 
tionless and lie at known and controllable distances from one 
another lfH-l5i l2lll . An alternative way is to identically prepare 
the atoms with no information about the orientation of their 
dipole moments. In this case, U12 should be averaged over 
all orientations of the dipole moments, so that cos 2 rj = 1/3 
and then 0x2 = at rx2 = A/4. The fact that O12 can be 
zero at rx2 = A/4 would allow for an experimentally difficult 
manipulation of atoms trapped at small distances. 

We first perform numerical analysis of entangled proper- 
ties of the system by solving the master equation (fTJ for the 
steady-state, p = 0, using the quantum optics toolbox for mat- 
lab 0711 . We quantify entanglement in terms of the concur- 
rence that relates entangled properties to the coherence prop- 
erties of the atoms pHl . 

Figure Q] shows numerical results for the concurrence plot- 
ted as a function of the Rabi frequency Oo for the case of 
resonant driving, = and different Ao and also different 
ratios of the atomic relaxation rates. We see that the concur- 
rence is insensitive to Oo until Oo = Ao- At this specific 
value of the Rabi frequency, the concurrence exhibits a sharp 
peak whose amplitude is independent of Aq. Consequently, 



the amount to which the atoms are entangled by dissipation 
is independent of the frequency difference Ao. Another in- 
teresting, and perhaps surprising feature seen in Fig. Q] is the 
dependence of the amplitude of the concurrence on the ratio 
72/71, that the seemingly symmetric system in respect to the 
rate of exchanging an excitation, i.e. 7x2 = 721, exhibits an 
asymmetry in the dependence of the concurrence on 72/71- 
The atoms are more entangled when the spontaneous emis- 
sion rate of the driven atom is larger than the other atom. 
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FIG. 1 : The steady-state concurrence versus the normalized Rabi frequency 
for Aj = (Q = flo) and different Ao and 72/71: Ao = 1571,72/71 = 
1 (solid line), Ao = 1571,72/71 = 5 (dashed line), Ao = 2571,72/71 = 
1 (long dashed line) and Ao = 2571 ,72/71 =5 (dash-doted line). 

We also see the concurrence plateau for Oo < Ao whose 
magnitude is independent of Oo but decreasing with increas- 
ing Ao- The presence of the plateau for small Ao indicates 
that at these frequency differences the atoms are indistinguish- 
able. This reflects the fact that photons emitted spontaneously 
by each atom induce in the other atom oscillations that are 
partially coherent with their own spontaneous oscillations and 
can thus lead to an entanglement. Thus, in order to entangle 
the atoms on demand, the frequency difference Ao should be 
quite large to prevent the atoms to exchange photons of an ar- 
bitrary frequency. In other words, the process of turning on or 
turning off entanglement on demand can be better controlled 
at large Ao. 

It is worth noting that similar behave of the concurrence is 
found for a detuned driving, / that the resonant peak 
appears at = A, where = (O^ + A|) 1/2 is the Rabi 
frequency in the detuned field and A = Ao + A^ . 

The numerical results have shown that the atoms decay to 
a stationary entangled state only when the Rabi frequency is 
tuned to the frequency difference between the atomic transi- 
tion frequencies, i.e. when O = A. We now calculate analyt- 
ically the density matrix of the system and the concurrence to 
explain the physical origin of this feature and to explore the 
importance of dissipation in the creation of the entanglement. 
In these calculations, we first diagonalize the Hamiltonian Hq 
of the system and find that it has eigenstates 



|1) 
|3) 



|ei>|- 



|2> = |ex}|- 
|4) = |.gi)|- 



(6) 
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with the corresponding energies 

E lt2 = ^h(A±n), £ 3 ,4 = ~ft(A=Ffi), (7) 

where 

|+) = cos0|e 2 )+sin%2}, h> = sin0|e 2 )-cos% 2 ), (8) 

are the dressed states of the driven atom and the angle 9 is 
determined by cos 2 9 = \ + Note that the energy 

states are in the form of product states indicating that the sub- 
systems are separable. Thus, the subsystems remain separa- 
ble until spontaneous emission is included. Furthermore, the 
states are non-degenerate unless O = A and then the states \2) 
and 1 3) become degenerate. In the physical terms, the condi- 
tion = A corresponds to resonance of the Rabi sideband 
ujl + O of the dressed-atom with the transition frequency of 
the undriven atom. 

We now include spontaneous emission and rewrite the 
master equation in terms of dressed-atom operators R~ = 
\-)(+\,R + = l+X-l. and R, = Kl+X+I - l-X-l), and 
make the unitary transformation on the density matrix of the 
system 



exp(iHt/h)p cxp(—iHt/h), 



(9) 



where H = hASf + KVtR z . 

The master equation for the transformed density operator is 
found to contain terms oscillating at frequencies A and A±fi. 
When fl 7^ A and the frequency difference is much larger 
than the atomic decay rates, A S> 71, 72, these terms oscillate 
rapidly in time and, in the limit of t — > 00, make negligible 
contribution to the master equation. After discarding them, 
which is a form of the secular approximation, we readily find 
that the master equation reduces to 



C\p = 71 (2S 1 p~S^~ — S^S-l p — p~S^~ S 1 



(10) 



which contains only the term representing the damping of the 
undriven atom. No terms representing the other atom are in- 
volved, which explains why entanglement is not found in the 
system when 51 7^ A. 

The situation changes when ft = A. In this case, the terms 
oscillating at A + CI become non-oscillatory and the other 
become rapidly oscillating at frequencies A and 2A. When 
we discard these rapidly oscillating terms, we find the master 
equation has the form 



p = C(p) = Cip + C d p + C c p, 



(11) 



where Cd is an operator representing the damping of the 
dressed-atom system 

£ d p = lQ ([Rz,pR z ] + [R z p,R z ])+j+([R-,pR + ] 

+ [R-p, R+] ) + 7 _ ( [R+, pR-\ + [R+p, R-]) , (12) 

with 70 = (72/4) sin 2 (29), 7+ = 72cos 4 0,7_ = 72 sin 9, 
and C c is the operator representing the dissipative coupling 
between the undriven atom and the dressed-atom system 



C c p = 712 ([pS+,R-] + [St, R-p] +H.c. 



with 712 = 712 cos 2 9. 

Here, the master equation contains the dissipative terms of 
both subsystems. However, no coherent terms are involved, 
indicating the absence of coherent dynamics in the system. 
Nevertheless, the system can be driven to a stationary entan- 
gled state due to the presence of the dissipative coupling be- 
tween the subsystems. 

In order to examine the final state of the subsystems, we 
solve the equation C(p) = to find the steady-state density 
matrix of the system. It is easily verified that in the density 
matrix, written in the basis of the product states (O, has a X- 
state form H 



(14) 



pn 
P22 P23 

P32 P33 
p 4 A 

where the non-zero density matrix elements are 

Pu = 7-7i 2 2 /A P22 = 7- (7i + 7+) 712/A 

P44 = 1 - 7- [77i (7 - 7o) + 37-712] / D > 

P23 = P32 =7i 7- (7-7o)7i2/A (I 5 ) 



with 



D = (7-7o) 2 [71 (7+ +7-) 



7i 2 2 ] 



+ (7+ + 7-) [7o7i (7 ~ 7o) + 47_7i2 



(16) 



and 7 = 71 + 70 + 7+ + 7- is the total damping rate. 

The arrangement that the laser field only couples to one 
of the two atoms resembles very much an open cascade sys- 
tem iRoL PTill . That is, the output (fluorescence) of the laser 
driven atom drives the other atom with no feeding back onto 
the first atom. In this case, the spontaneous emission coupling 
term has the form 



£ c p = 7i2 {[St,R-p] + [pR + ,Si]) , 



(17) 



and then it is straightforward to show that the steady-state den- 
sity matrix elements are very similar in form to the ones found 
for the mutually driving systems. They differ only in that the 
denominator D is replaced by 



D' = (7+ + 7_) [771 (7 - 7o) + 27-7 2 2 



(18) 



(13) 



The X-state form of the density matrix allows us to obtain a 
simple formula for the concurrence quantifying the amount of 
entanglement between the undriven and the dressed-atom sys- 
tems at the level crossing point of SI = A. The concurrence is 
shown in Fig.|2]for continuously varying 72/71 and A^. The 
left frame illustrates the concurrence for the case when both 
subsystems react to the mutually emitted photons, whereas the 
left frame illustrates the concurrence for the case of an open 
cascade system, where the undriven atom reacts to the photons 
emitted by the dressed-atom system, while there is no interac- 
tion in the reverse direction. From the data, we note that the 
concurrence reaches a maximum height for a large positive 
laser detuning, Al > 0, and 72 > 71. In terms of the pop- 
ulation of the dressed-states of the driven atom, the condition 
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Al > corresponds to no inversion on the dressed-atom tran- 
sition, the higher frequency Rabi sideband, which for SI = A 
is resonant to the transition frequency of the undriven atom. 
The atoms remain separable for practically all 71 > 72. Thus, 
the driven atom must decay faster than the undriven one to en- 
tangle the atoms by dissipation. Needless to say, the fluores- 
cence from the driven atom acts as a pump for entanglement 
between the atoms even if both systems react to the mutually 
emitted photons. 




FIG. 2: Steady-state concurrence as a function of a = (71 — 72 )/ (71 + 72 ) 
and cos 2 9 for mutually driving subsystems (left frame) and for the case of a 
cascade open system (right frame). 

The shape of the concurrence for the cascade open system, 
illustrated in the right frame of Fig. [2] continues the trend of 
the case of the mutually driving systems, shown in the left 
frame. However, it should be pointed out that the analogy 
between the two cases is not absolute because there is always a 
non-zero probability that the driven atom will absorb a photon 
emitted by the undriven atom. Despite this, Fig. [2] shows that 
the concurrence is formally almost the same for both cases 
and therefore the system can be regarded as an open cascade 
system. 



A non-zero value of the concurrence reveals the existence 
of entangled states in the system. These can be easily found 
diagonalizing the density matrix (TBI) , which results in super- 
position (entangled) states 

|s) = cos0|2) +sin0[3), \a) = sin0|2) -cos0|3), (19) 

where the angle cf> is determined by cos 2 <f> = ^+<V with 

S = P22—P33 an d G = \/ 5 2 + 4|p23| 2 - Clearly, the stationary 
state of the system is a mixed state involving two entangled 
states. The mechanism for the creation of the entangled states 
is the coherence /?23 generated by the interaction of the atoms 
with a dissipative reservoir. 

In summary, we have proposed a scheme for the creation of 
entanglement by the dissipative process of spontaneous emis- 
sion. The scheme involves two non-identical two-level atoms 
separated at a small distance and interacting with a common 
reservoir. Arranging the system such that the undesirable 
dipole-dipole interaction between the atoms is suppressed and 
only one atom is addressed by an external laser field, we have 
found that the atoms can be entangled on demand by tuning 
the Rabi frequency of the laser field to the difference between 
the atomic transition frequencies. At this frequency, two en- 
ergy levels of the system become degenerate resulting in a 
large coherence causing the system to evolve into a stationary 
entangled state. We have also found that the entanglement oc- 
curs only when the damping rate of the driven atom is larger 
than the other atom. We have demonstrated that this system 
predominately behaves as a cascade open system where the 
fluorescence from the dressed atom drives the other atom with 
no interaction in the reverse direction. 
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